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ABSTRACT 


The  rapid  advancement  of  technology  in  the  electronic 
industry  has  created  a  need  for  solving  a  certain  class  of 
large  nonlinear  reliability  problems.  The  Deportment  of 
Defense  and  civilian  industries  when  planning  for  the 
allocations  of  budget  resources  require  an  efficient  method 
for  solving  these  problems.  In  certain  areas  a  specific 
reliability,  such  as  with  nuclear  reactors,  must  be 
attainable  to  ensure  the  accomplishment  of  the  mission  or 
to  provide  necessary  safety  factors.  The  current  methods 
used  to  solve  these  problems  are  very  difficult  and  time 
consuming.  A  computer  algorithm  which  is  based  on 
geometric  programming  has  been  developed  to  provide  a 
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Chapter  1 
INTRODUCTION 


Problem  Statement 

In  the  general  field  of  electronic  instrumentation,  an 
important  goal  for  both  government  and  private  industry  is 
that  of  fielding  systems  that  can  operate  reliably  for  a 
specified  period  of  time  within  specified  environments. 

One  of  the  most  difficult  problems  in  this  field  is 
ensuring  a  level  of  acceptable  reliability  at  minimum 
cost.  Related  to  this  problem  is  its  economic  dual,  which 
is  to  maximize  reliability  within  fixed  cost  limits.  The 
purpose  of  this  thesis  is  to  show  that  these  types  of 
reliability  problems  are  easily  solvable,  within  certain 
restrictions,  using  geometric  programming  (GP). 

Component  reliability  is  defined  as  the  mathematical 
probability  that  an  individual  component  will  continue 
performing  its  Intended  mission  for  a  specific  length  of 
time.  The  failure  rate  of  a  particular  component  is  stated 
as  a  probability  between  0  and  1.  In  mathematical  terms, 
component  reliability  equals  1  minus  probabi 1 ity  of 
failure,  or  1  -  P( failure).  If  a  system  is  composed  of 
more  than  one  component  in  series,  then  the  system 


reliability  would  be  the  product  of  the  individual 


component  reliabilities.  For  one  type  of  component,  the 
system  reliability  would  become 
1  -  P(fallure)*  (of  un,ts>. 

Although  the  definition  of  component  reliability  is 
straightforward,  systems  involving  reliability  can  become 
rather  complex.  Figure  1.1  depicts  a  general  reliability 
system.  The  system  consists  of  n  elements  where  each 
element  can  have  one  or  more  redundant  or  parallel 
components. 

The  reliability  of  any  separate  element  and  the 
reliability  of  the  system  itself  are  increased 
exponentially  through  redundancy  of  components.  However, 
one  must  then  face  the  correspond i ng  mu  1 1  f  p 1 1 cat i ve 
i ncrease  i n  cost . 

The  first  problem  of  achieving  acceptable  reliability 
at  minimum  cost  could  be  stated  as  follows.  Consider  a 
system  involving  electronics  with  n  interconnected 
elements.  Each  element  Is  made  up  of  some  number  of  backup 
components  that  will  switch  on  in  case  of  failure.  Figure 
1.1  is  such  a  system.  A  two-element  system,  where  we  wish 
to  minimize  the  cost  of  paral lei  back-up  components  subject 
to  a  minimally  acceptable  level  of  total  system 


reliability,  would  be  formulated  as 
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minimize  cost:  C  =  $5Nj  +  $7N2 
subject  to:  \l  -  0.4  1 


>  0.90  (1.1) 


where  Nj  represents  the  number  of  components  of  element  1 
and  represents  the  number  of  components  of  element  2. 

The  cost  of  each  Nj  component  is  $5;  for  each  component 
it  is  $7.  The  probability  of  failure  of  each  component  of 
element  1  is  0.4;  the  probability  of  failure  of  each 
component  of  element  2  is  0.6;  while  the  minimally 
acceptable  level  of  total  system  reliability  is  0.90. 

The  second  problem,  the  economic  dual,  also  considers  a 
system  involving  electronics  with  n  interconnected 
elements.  As  before,  each  element  is  made  up  of  a  number  of 
backup  components  that  will  switch  on  in  case  of  failure. 
However,  in  the  case  of  the  two-element  system  shown  below, 
the  problem  is  maximization  of  total  reliability  given  a 
total  cost  constraint,  e.g.,  $60. 


maximize  reliability:  R 
subject  to: 


1 


N 

0.4 


N. 


0.6 


$5Nj  +  $7N2  <  $60  (1.2) 


This  thesis  will  develop  an  algorithm  for  the  solution 
of  both  types  of  problems.  The  algorithm,  using  geometric 


r: 


it- 


& 
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programming,  will  be  programmed  for  interactive  use  in 
Basic,  and  a  class  of  appropriate  test  problems  will  be  used 
to  demonstrate  the  algorithm's  ability  to  achieve  global 
optimal ity. 


The  General  Geometric  Programming  Approach 

Prob 1 ems  i nvo 1 v i ng  reliability  are  norma  1 1 y  charac¬ 
terized  as  nonlinear  design  problems.  The  objective 
function  and/or  the  constraint  functions  take  the 
genera 1 i zed  non 1 i near  form  of 

minimize  g<X)  =  S  C  P  (X) 

i  =  l  i  i 


where 

3 1  |  ^  j  „  3* 

P.(X)  *  X,  X_  ...  X  1  =  I,  2 . n; 

1  i  c  n 

Cj  and  a^j  are  typically  physical  constants;  and 
Xj  are  design  variables,  J  =  l,  2,  ....  n. 

Geometric  programming  (GP)  is  a  mathematical 
programming  technique  used  to  solve  such  nonlinear  design 
problems.  The  general  form  of  a  GP  problem  is 

minimize  cost:  C  =  gQ(X) 


-'v* 


."■'.y.'-.v."-' 
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subject  to: 


where 


g, (X)  <  1 


XJ  >  0 


f  1*  2 ,  .  .  .  ,  fn 


j  —  1*  2 .  •  •  •  *  n 


T  a 

I  n  I  jt 

g  (X)  *  r  C  r  X 

1  t= 1  ft  j=  1  j 


A  few  definitions  used  in  GP  are  helpful  in  order  to 
further  understand  and  discuss  the  algorithm.  A  "term"  is 
defined  as  a  grouping  of  one  or  more  variables  with  or 


without  coefficients  separated  from  any  other  grouping  of 
variables  by  a  plus,  minus,  or  inequality  sign.  The  degree 
of  difficulty  is  an  i nd i cat l on  of  how  difficult  it  is  to 
solve  the  original  problem;  progressively  higher  numbered 
problems  become  more  difficult  to  solve.  The  degree  of 
difficulty  (DO)  of  a  particular  problem  is  determined  by 
taking  the  number  of  individual  terms  in  the  entire  problem 
minus  the  number  of  variables  minus  1.  Thus,  a  GP  problem 
such  as 


minimize  cost: 


C  =  *12N  9N_ 

i _ ! — £j 

term  i 


subject  to: 


Nj  >0.95 
term  2 
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N  >  0.95 

L  — _ i 

term  3 


would  have  a  DO  of  zero,  as  3  terms  minus  2  variables 
(Nj  and  N^)  minus  l  equals  zero. 


In  GP,  a  zero  DD  problem  is  desirable,  as  a  higher  DD 
problem  would  require  some  type  of  mathematical  manipulation 
before  the  problem  could  be  solved.  The  reliability 


problems  solved  using  the  algorithm  developed  below  will  be 
one  DD,  and  therefore  will  require  an  internal  mathematical 
method  to  solve  for  the  unknown  6  values.  This  method  will 
be  addressed  in  Chapter  3.  The  first  step  in  developing  the 
required  algorithm  is  to  rewrite  the  problem  using  various 
substitutions  to  manipulate  the  variables  so  that  the 


resulting  GP  problem  is  one  DD. 


.V.  «\V,  c,  .*.  /  /  V  ■/  ,  . 
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Chapter  2 

CURRENT  METHODS  OF  SOLVING 
RELIABILITY  PROBLEMS 


8 


Reliability  as  a  separate  mathematical  science  was 
developed  as  a  direct  result  of  analyses  conducted  during 
World  War  II.  The  rapid  automation  of  technology  created  a 
need  for  a  way  to  determine  the  reliability  of  specific 
items  of  equipment  as  well  as  total  systems. 

At  the  close  of  World  War  II.  the  military  conducted 
several  studies  on  reliability.  These  studies  revealed 
some  startling  facts:  (1)  electronic  equipment  used 
during  Navy  maneuvers  was  operational  only  30X  of  the  time; 
(2)  67-75 X  of  the  Army  electronic  equipment  was  out  of 

commission  or  under  repairs;  (3)  Air  Force  repair  and 
maintenance  costs  were  ten  times  the  original  equipment 
costs;  (4)  for  every  electonic  tube  on  the  shelf,  seven 
were  in  transit;  (5)  approximately  one  electronic 
technician  was  needed  for  every  250  electronic  tubes:  and 
(6)  the  number  of  electronic  tubes  required  on  board  a 
destroyer  had  risen  from  60  to  3,200  (Shooman,  I960). 

In  1950,  the  Department  of  Defense  established  an  ad 
hoc  committee  on  reliability,  which  was  replaced  In  1952  by 


the  Advisory  Group  on  the  Reliability  of  Electronic 
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Equipment  (AGREE).  AGREE  and  the  Institute  of  Electrical 
and  Electronic  Engineers  (IEEE)  are  currently  the  main 
agencies  conducting  reliability  analysis. 

Today  reliability  is  stressed  in  military  applications 
of  electronics,  aviation,  and  weapons  systems.  Reliability 
has  become  an  important  aspect  of  military  contracting  as 
it  is  essential  in  evaluating  the  system  usefulness  or 
goodness  along  with  cost,  size,  and  weight. 

The  approach  to  solving  reliability  problems  varies 
with  each  individual  and  organization.  Most  often, 
reliability  problems  are  treated  as  applied  probability  or 
statistics  problems.  Because  analysis  of  complicated 
reliability  problems  becomes  lengthy  or  difficult, 
transform  methods  or  computer  solutions  are  necessary. 
Complex  systems  are  almost  impossible  to  solve  in  entirety 
and  must  be  decomposed  into  functional  areas. 

Currently,  reliability  problems  of  the  type  described 
in  Chapter  1  are  theoretically  solved  using  the  Kuhn-Tucker 
conditions.  In  reality,  most  reliability  problems  are 
solved  using  either  a  variable  search  or  Markov  process 
method.  It  is  fairly  simple  to  solve  reliability  problems 
using  the  Markov  process;  however,  the  labor  involved 


increases  tremendously  as  the  number  of  elements  in  the 
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total  system  expands.  In  addition,  this  method  cannot 
guarantee  global  optimality;  geometric  programming,  on  the 
other  hand,  does  guarantee  global  optimality. 

GP  is  practical  for  solving  reliability  problems  since 
it  is  capable  of  exploiting  the  linear  algebraic  structure 
of  each  problem.  These  linearities  can  appear  as  linear 
equations,  linear  inequalities,  or  as  matrices  associated 
with  non  1 i near i ties. 

Research  on  the  type  of  reliability  problems  discussed 
in  Chapter  1  revealed  only  one  previous  attempt  to  use  GP. 
In  1968,  A.  J.  Federowi cz  and  M.  Mazumdar,  while  working 
for  the  Westinghouse  Research  Laboratories,  developed  a 
method  of  using  substitution  to  enable  the  given  problem  to 
be  written  in  standard  GP  form.  Once  the  problem  was 
written  in  GP  form  it  was  easier  to  solve. 

The  algorithm  developed  in  this  thesis  follows  the 
initial  substitution  and  manipulation  that  Federowi cz  and 
Mazumdar  used  to  solve  maximization  problems.  Instead  of 
using  Federowi cz  and  Mazumdar 's  method,  which  uses 
logarithms  and  di fferent fat fon  to  find  the  unknown  ; 
values,  either  the  quadratic  equation  or  the  Brent-Oekker 
method  was  used  here  to  find  the  positive  root  of  the 
unkown  i.  A  method  for  solving  minimization  problems  has 
a  I  so  been 
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Included.  Further  comparison  between  Federowicz  and 
Mazumdar's  method  and  the  one  developed  in  this  thesis  will 
be  addressed  in  Chapter  5. 
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Chapter  3 
METHOO  FOR  SOLVING 
MINIMUM  COST  PROBLEMS 

There  are  four  phases  In  the  GP  method  for  solving 
minimum  cost  problems  and  in  developing  an  associated 
generalized  computer  algorithm.  Manipulation  of  terms  and 
substitution  in  the  early  phases  can  reformulate  the 
original  problem  in  the  standard  GP  form  used  by  Woolsey's 
"Quick  &  Dirty"  four  rule  algorithm  (see  Appendix  A). 

Phase  1:  The  Objective  Function 

This  phase  consists  of  three  steps  and  uses 
substitution  and  the  rules  of  logarithms  to  manipulate  the 
objective  function  of  the  original  problem.  In  the  case  of 
the  first  problem  discussed,  minimizing  cost  subject  to  a 
specified  reliability,  the  generalized  form  of  the 
objective  function  is 

n 

minimize  cost:  C  =  S  C  N 

j  =  l  J  J 

where 

j  *  I,  2,  ...,  n 

To  transform  the  objective  function  into  the  form 
required  by  GP .  consider  the  following  procedure. 
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Step  1)  Using  substitution,  let 


N.  «  in  X. 


Step  2)  Rewrite  the  objective  function  as 


n 


minimize  cost:  C  =  £  C  In  X 

j=l  J  J 


Step  3)  Use  the  rules  of  logarithms  to  get  the 
objective  function  in  standard  GP  form. 

a)  When  the  anti  log  of  a  function  such  as 

in  Yj  =  Cj  in  Xj 

i s  taken ,  i t  becomes 


b)  When  the  anti  log  of  an  additive  function 
such  as 


In  Y  =  C j  In  X j  +C2  In  X? 


is  taken,  it  becomes  multiplicative 


Y 


Transform  the  original  objective  function  of  sample 
problem  (1.1)  by  performing  the  above  steps. 


;  .'  T 


«  -  *  *  .  *  «rM  *  ,  i  ,  •  - 

*w*.  .v.  a*./-* 
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Step  0)  Restate  the  original  objective  function, 
minimize  cost:  C  =  $5Nj  +  $7^ 

Steps  1  &  2)  Rewrite  the  above  equation  as  discussed 

minimize  cost:  C  =  5  In  Xj  +  7  In 

Step  3)  Take  the  anti  log  of  the  above  equation. 

C  5  7 

minimize  cost:  e  - 

The  objective  function  is  now  in  standard  GP  form. 
Phase  2:  The  Constraints 

This  phase  consists  of  nine  steps  and  uses  substi¬ 
tution  and  the  rules  of  logarithms  to  manipulate  the 
constraint  of  the  original  problem.  The  general ized  form 
of  the  orignal  constraint  discussed  is 


Step  1)  The  reliability  o'  the  jth  element  stated  as 


1 


■ aj 
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is  difficult  to  handle  in  GP.  Use  substitution  to 
get  rid  of  this  "messv"  term  bv  letting 


Z 


j 


N 


j 


Step  2)  Rewrite  the  constraint  as 


z,Z2-.-zn  >  b 


SteD  3)  At  this  ooint,  because  of  the  substitution 
oreviously  performed,  the  objective  function  is 
written  in  terms  of  X.,  while  the  constraint  is 
written  in  terms  of  Z..  The  problem  cannot  be  solved 
in  this  form,  as  there  is  now  no  relation  between  the 
variables  in  the  objective  function  and  the 
constraint.  In  GP.  all  constraints  will  be  tight  at 
optimality  (Duff  in  et  a  I..  1967).  Therefore,  an 
additional  constraint  (written  as  an  ineaualitvl  can 
be  added  for  each  Z..  showing  the  substitution  used  in 
Step  1  . 

NJ 

Z  .  <  1  -  P  .  ] 

i  —  l 


Step  4)  Terms  that  have  a  constant  raised  to  a 
variable  power  are  rather  difficult  to  handle  with 
GP.  Thus,  let 
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N  . 

Y  .  =  P  .  • 

I  1 


?J| 


Step  5)  Take  the  natural  log  of  both  sides 


In  Y  =  N  In  P 


Step  6*1  Use  the  substitution  N.  =  In  X.,  as  was  done 
for  the  objective  function,  to  obtain 


I n  Y  .  =  In  X  .  In  P . 


or 


In  Y  =  In  P  j  In  X  ( 


Step  7)  Take  the  anti  log  of  the  above  eauation. 


JF1 

TS 

•;v 


W  1 

% 
*  ' 

£ 


In  P 


Y1  "  Xi 


Step  8)  Rewrite  each  added  constraint  (as  reauired 
in  Step  3).  using  this  manipulation  to  obtain 

In  P  . 

Zj  <  .  -  X.  1 

Step  9)  Manipulate  each  constraint,  so  that  the 
riaht-hand  sides  are  less  than  or  eaual  to  l. 


a)  Dividina  both  sides  bv  Z 


constraint  becomes 


i  >  b| Z. 1 Z" 1 . . .Z'1 I 
1  Z  n  i 


the  first 
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b)  Rewritten.  It  becomes 


I 


c)  The  added  constraints  become 

In  P, 


Zt  +  X 


J 


<  1 


By  performing  the  above  nine  steps  to  the  constraint 
of  the  original  sample  problem  (1.1),  the  resulting 
constraints  are  transformed  into  standard  GP  form.  An 
application  of  the  above  nine  steps  is  shown  below. 


Step  0) 


Restate  the  original  constraint: 


>  0.90 


Steps  1  &  2)  Replace  the  "messy"  terms  with  Zy 
and  rewrite  the  constraint: 


ZjZ2  >  0.90 

Steps  3-8)  Add  the  additional  contra ints  for  each  Z 
as  discussed: 


Z 


2 


> 


In  0.6 
2 
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Steo  9)  Rewrite  each  of  the  constraints  as  discussed. 


O.SZj lZ2l  <  I 


z,  +  x["  °*4  <  1 


Z2  +  *2°  °'6  i  1 


The  constraints  are  now  written  in  standard  GP  form. 


Phase  3:  Solving  bv  GP 

This  Dhase  consists  of  seven  steps,  and  uses  Woolsev's 
four  rules  to  solve  the  problem.  Phases  1  and  2  have 
reformulated  sample  problem  (1.1)  as  follows. 


minimize  cost:  C  =  X^  X^ 

term  1 


subiect  to: 


0.90Z  [Z~X  <  1 

* _ ~J 

term  2 


Z,  *  *!"  °-4  £  I 
term  3  term  4 


,  x^"  °-6  <;  l 


term  5  term  6 


At  this  point  the  resulting  GP  problem  is  one  DO  as 


there  are  6  terms  and  4  variables  ( X j ,  X^.  Zj ,  and  Z^). 
Using  Woolsey's  tour  rules  ot  geometric  programming 
(Appendix  A),  we  can  solve  sample  problem  (1.1). 


Step  1)  Write  the  form  of  the  optimal  solution 
according  to  Rule  1. 


minimize  cost:  C 


(5 


3 


+ 


=  (1/5 l1 (0.9/:2)2(i2 
+  5  ) 

4  (  1  /  5  5 ) 5  (  1  /  *  6  )  6  (  5  5 


)2( I /: 3 ) 3 ( I / i 4 


*4 

) 


Step  2) 

State 

the  exponent 

matrix  of  the 

S's 

accord i 

ng  to  Rule  2. 

(OF)  : 

i 

'  1 

=  1 

(3.1) 

(Xj  )  : 

55 1 

+  InO. 4‘ . 

=  0 

(3.2) 

(X£)  : 

751 

+  1 nO . 6?6 

=  0 

(3.3) 

(Zj  )  : 

-*2 

+  '  3 

=  0 

(3.4) 

(Z2)  : 

_  £ 

2 

+  S5 

=  0 

(3.5) 

Step  3)  Solve  for  the  values  of  the  $'s  using  the 
logic  that  follows.  A  pattern  occurs  in  the  i  matrix 
table.  Regardless  of  the  number  of  variables  in  the 
original  problem,  this  pattern  will  always  be 
present.  The  pattern  is  shown  below. 


vs 
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f)  Thus 


;  _  i.  s  ‘ 

2  '  3  5 


Step  4)  Skfp  to  Rule  4  since  the  values  of  >. 
and  are  unknown. 


a)  Write  the  equations  for  *  and  *5- 


s3  =  zl(f3  +  s4) 


*5  -  Z2(i5  +  ' 6> 


b)  Solve  each  equation  in  terms  of  Z 


Z1  “ 


(53) 


(<3  +  *4> 


Z2  - 


(  ‘  5  > 


(f5  +  V 


c)  Replace  the  S's  with  the  known  values 


21  " 


(  ^2) 


(«  +  5.4567833) 


Z2  " 


<  '  2 ) 


(  %  +  13.703306) 


-::v 


V_’.-  V. V-’ 
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b)  Substitute  in  the  values  found  in  SteD  4e 
and  solve  for 

-(-17.244081+  [(  17.24408)2-4(  .  1  )  ( -67 . 29837  )] 

2  2(  .  n 

Thus , 


*  =  176.25895 

and 

i 2  =  -3.818155 

c)  Since  In  GP  the  :'s  must  be  nonnegative,  onlv 
the  rea I  pos i t i ve  root  of  the  quadrat i c  eguat i on 
can  be  used;  therefore,  drop  the  negative  : 

i2  =  :3  =  ?5  =  176.25895 

d)  If  the  power  of  the  first  term  is  greater 
than  2.  then  to  solve  for  the  unknown  *  a 
modified  secant  method  may  be  used.  This 
method,  known  as  the  Brent-Dekker  method,  is 
based  upon  both  bisection  and  the  secant  rule 
methods.  It  is  similiar  to  the  Newton-Raphson 
method,  except  only  the  ability  to  evaluate  f(x) 
is  necessary.  The  Newton-Raphson  method  is 
slightly  faster,  but  it  reguires  the  first 


■  •  Vs 


v  ‘  .vv  v  v  •>  ■ . v 


1  *jf  mJCm  m,  .1* kjCajCkjC. 4V«V. 


*  ■*  fif.Ajf. 


approximation  concerning  the  bracketing  of  the 


roots  to  be  close.  In  addition,  the  Newton- 
Raphson  method  may  not  converge,  whereas  the 
Brent-Oekker  method  i s  guaranteed  to  work  once 
the  function  is  bracketed. 

To  start  the  Brent-Dekker  method,  attempt 
to  bracket  the  function  with  two  values  (B  and  C) 
Figure  3.1  depicts  this  starting  bracket.  The 
function  must  be  continuous,  and  f(B)f(C)  must 
be  less  than  zero.  Throughout  this  algorithm,  B 
Is  assumed  to  be  the  better  root. 

Convergence  is  based  on  a  mixed  relative- 
absolute  error  test  where 

C  -  B 

-  £  max[ Absolute  Error, tfi1*  Relative  Error] 

2 

The  secant  rule  calculates  the  next  iterate 
(D)  by  starting  with  the  first  two  iterates  A  and 
B .  The  var table  A  is  initially  set  equa 1  to  C . 
Thus , 

B  -  A 

0  =  B  -  f (B )  - 

f ( B)  -  f ( A ) 

Figure  3.1  depicts  how  the  starting  bracket  is 
used  to  find  the  next  iterate  D  using  this 


method.  If  0  is  to  become  the  next  iterate,  thei 
it  must  lie  between  B  and  the  midpoint  between  B 
and  C.  If  0  is  between  B  and  the  midpoint,  then 
it  is  set  equal  to  B.  The  old  B  is  then  set 
equal  to  C.  If  D  is  not  in  the  interval  between 
B  and  the  midpoint,  then  the  midpoint  between  B 
and  C  is  used  as  the  new  iterate. 

The  Brent-Dekker  method  guarantees  that 
there  is  either  a  root  of  f(x)  =  0  in  the 
Interval  between  B  and  C,  or  that  one  of  the  end 
points  is  as  close  to  a  root  as  the  desired 
precision  permits,  if  the  function  being 
evaluated  is  continuous  and  the  original  bracket 
(B,C)  results  in  f(B)f(C)  <  0. 

Step  6)  Solve  for  Xj  and  according  to  Rule  4,  by 

replacing  the  5's  with  their  known  values. 

a)  Write  the  equations  for  S.  and 
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b)  Substitute  in  the  known  values  for  the  ;'s. 


5.4567833  =  x  0 * 9 1 62907 n 76 . 25095  +  5.4567833) 


13.703306  =  X2°‘5108256( 176.25895  +  13.703306) 


c)  Solve  the  above  equations  for  Xj  and  X2 . 


Xj  =  45.871487 


X2  =  171.904510 


Step  7)  Use  the  above  values  to  solve  for  the 
ob j ect i ve  f unct i on . 


a)  Find  the  optimal  GP  values  of  and  N2  by 
substituting  in  the  values  for  the  In  Xj. 


Nj  =  In  Xj  =  In  45.871487  =  3.8258437 


N2  =  In  X2  =  In  171.90451  =  5.1469392 


b)  Therefore,  the  optimal  GP  value  of  the 


objective  function  is 


minimum  cost:  C  =  $55.16 


,v*'  "v  *►  •%' 
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Phase  4:  Integer  Solution 

Components  cannot  be  split  into  fractions,  therefore 
an  integer  solution  is  required.  This  phase  will  consist 
of  four  steps  that  will  ultimately  provide  the  minimum  cost 
integer  component  solution. 


Step  1)  Round  each  Nj  down  to  the  nearest  integer, 
then  add  1.  and  write  its  value. 


N.  =  4 


N2  =  6 


Step  2)  This  value  will  more  than  satisfy  the 
original  constraint;  however,  it  may  not  be  the 
optimal  combination  of  variables  that  both  satisfies 
the  constraint  and  has  the  minimum  cost.  Determine 
the  reliability  and  cost. 


reliability  =  1  -  0.44  1  -  0.66  =  0.9289384  >  0.90 


cost  =  $5(4)  +  $7(6)  =  $62.00 


Step  3)  Determine  the  number  of  combinations  that 
must  be  compared  to  ensure  that  the  optimal  solution 
i s  found . 


A  V-  A  A  .  r*.  A.  A  A.'  A  A  A  A  A  -  .  A  -  V.  A  •  •*“*-*'*- 


a)  Add  1  to  each  integer  value  found  in  Step  1. 
Determine  the  new  reliability  for  each  element. 


reliability 

element  1  = 

1  -  0.45 

=  0.989760 

reliability 

e 1 ement  2  = 

1  -  0.67 

=  0.9720064 

b)  Use  the  following  formula  developed  by 
Steven  A.  Strauss  to  find  the  number  of 
combinations  to  be  checked  to  ensure  that  the 
minimum  cost  is  found  subject  to  the  minimum 
required  reliability. 

+  combination  Nj  =  +  Integer(Nj)  -  integer 

C 1 n C 1  -  minimum  required  rel iabi 1 ity/k) ] 

of  -  -  — - - - - - - 

In  [ P ( f a i 1 ure  of  j ) ] 

where 


k  =  reliability  of  each  element  found  in 
Step  3a,  times  all  the  other  element 
reliabilities  except  for  the  jth  element 
(the  one  whose  optimal  number  of  combin¬ 
ations  is  currently  being  calculated) 


Therefore 
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+  combination  Nj  =  +  4  -  integer 


[in  (1  -  0.90/0.9720064)] 
of  -  =  +  2 


In  0.40 


±  combination  =  ±  6  -  integer 


[in  (1  -  0.90/0.9897600)] 

of  -  =  +  2 

in  0.60 

c)  Find  ail  combinations  for  each  variable  as 
determined  in  Step  3b.  The  combinations  are 
also  a  sensitivity  analysts  for  reliability  and 
cost.  Table  3.1  shows  the  combinations  of  N( 
and  which  result  for  example  problem  (1.1). 


Step  4)  Select  the  feasible  combination  with  the 
minimum  cost.  This  will  be  the  optimal  integer 
solution. 


minimum  cost:  C 


N 

N 


1 

2 


$60.00 

5 

5 
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Table  3.1  Minimum  Integer  Combinations 


N. 

N2 

RELIABILITY 

COST 

(IN  DOLLARS) 

2 

4 

0.7311360* 

38.00 

2 

5 

0.7746816. 

45.00 

2 

6 

0.8008090. 

52.00 

2 

7 

0.8164854. 

59.00 

2 

8 

0.8258912 

66.00 

3 

4 

0.8146944* 

43.00 

3 

5 

0.8632166. 

50.00 

3 

6 

0.8923300 

57.00 

3 

7 

0.9097980 

64.00 

3 

8 

0.9202788 

71.00 

4 

4 

0.8481 1 78* 

48.00 

4 

5 

0.8986307 

55.00 

4 

6 

0.9289384 

62.00 

4 

7 

0.9471230 

69.00 

4 

8 

0.9580338 

76.00 

5 

4 

0.8614871  * 

53.00 

5 

5 

0.9127963 

60.00 

5 

6 

0.9435818 

67.00 

5 

7 

0.9620531 

74.00 

5 

8 

0.9731358 

81.00 

6 

4 

0.8668348* 

58.00 

6 

5 

0.9184625 

65.00 

6 

6 

0.9494391 

72.00 

6 

7 

0.9680251 

79.00 

6 

8 

0.9791766 

86.00 

*  denotes  those 

combinations 

which  are  not 

feasible  since  they  are  less 
minimum  required  reliability. 

than  the 
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Chanter  4 


METHOD  FOR  SOLVING 


MAXIMUM  RELIABILITY  PROBLEMS 


There  are  three  phases  in  the  method  for  solving 
maximum  reii ability  problems  and  in  developing  an 
associated  generalized  computer  algorithm.  This  chapter 
uses  the  same  manipulation  of  terms  and  substitution  for 
the  objective  function  and  constraints  as  does  Chapter  3. 
The  second  problem,  maximizing  reliability  subject  to  a 
spec i f i ed  cost  1 i m i tat i on ,  wou Id  be  written  in  genera  I i zed 
form  as 


maximize  reliability:  R  = 


subiect  to: 


where 


N, \  n2 

1  -  P,  I  -  P/ 

*  C- 


:  c  n  <  b 
j=l  I  I  ~ 


j  =  1  .  2 . 4 


I  -  P 


Phase  1 :  Standard  GP  Form 

This  phase  consists  of  three  steps  used  to  get  the 
original  samole  problem  (1.2)  into  standard  GP  form.  This 
is  accomplished  bv  using  the  subst i tut  ions  and  manipu¬ 
lations  described  in  Chapter  3.  The  substitutions  include 
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N .  =  tn  X  . 
i  1 

and 

In  P. 

2J  =  1  ~  Xi  ’’ 

Step  1)  Restate  the  original  samole  problem  (1.2) 


f  Nl] 

1  nJ 

maximize  reliability:  R  =  ^1  -  0.4 

(  1  -  0.6  2 

subject  tot  SSN^  +f7N^  f_  $60 


Steo  2)  Use  the  substitution  and  manipulation  as 
described  in  Chanter  3  to  rewrite  samole  problem 
(1.2). 


maximize  reliabilitv:  R  =  ZjZ., 


subject  to: 


y  Sy  7  ,  60 

X  j  *2  —  e 


Z,  <  1  -  xjn  °*4 


z2  <  1  -  X-  0,6 


Steo  3)  GP  only  handles  minimization  problems: 
therefore,  to  change  the  obiective  function  to  a 
minimization,  invert  the  function.  The  resulting 
orob I em  becomes 


"•r.T'Af  A.  vv  ."JV"  «-.  ■ 

-"■Ar-a  AiV.  w.  /-  .--V-  .-j  ,Cvl>»-  .v  ^V- x- 
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minimize  reliabilitv:  R 


subiect  to: 


Z  1 Z  1 
1  2 

term  1 


term  2 


2  4  x'n  °-4  /  i 

i _ Lj  < _ ! _ i 

term  ?  term  a 


z2  +  x’n  0,6  <  I 

term  5  term  6 


Phase  2:  Solving  bv  GP 

This  ohase  consists  of  sW  steDs  and  uses  Wooisev's 
four  rules  to  solve  the  problem.  Sample  Droblem  (1.2)  is 
now  in  standard  GP  form  and  is  a  one  DO  oroblem  as  there 
are  6  terms  and  4  variables  (X } ,  X^.  Z( .  and  ) .  The 
following  steps  are  the  same  as  thosed  used  in  Phase  3  of 
Chapter  3 . 


Step  l)  Write  the  form  of  the  optimal  solution 
according  to  Rule  I: 
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5  5  ? 

minimize  reliability:  R  1  =  (l/ij)l(e  (l/e^)3 


i  ( S  +  i  )  5  5  ( i 

(1/-54)4(S3  +  -4)  3  4(  l/45)5(  l/i6)6(i5  +  S6)  5 


> 


Step  2)  State  the  exponent  matrix  of  the  5's 
according  to  Rule  2. 


(OF)  : 

:  1 

(Zj  )  : 

1  +  5  3 

(Z2)  : 

“SI  +  *5 

(Xj  )  : 

+  5&2  +  In  0 . 4 S 4 

(X2)  : 

+  752  + 

I  (4.1) 
0  (4.2) 

0  (4.3) 

0  (4.4) 

0  (4.5) 


Step  3)  Solve  for  the  values  of  the  i's  using  the 
logic  that  follows.  A  pattern  occurs  in  the  5  matrix. 
Regardless  of  the  number  of  variables  in  the  original 
problem,  this  pattern  will  always  be  present.  The 
pattern  is  shown  below. 


For  j  =  1  to  the  number  of  variables 
(2* J  +  1 )  1  1 

i(2*j  +  2)  =  (“Costj  /  In  Failurej)?2 

The  following  steps  depict  how  the  pattern  works  for 
the  two- variable  case. 


Iv-tvi-X-iv 
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a)  From  equation  (4.1)  determine  that 


b)  From  equation  (4.2)  determine  that 


c)  From  equation  (4.3)  determine  that 

'  5  -  1 

d)  From  equation  (4.4)  determine  that 

S2  =  0. 1832581 i4 

e)  From  equation  (4.5)  determine  that 

5  2  =  0.0729751  56 

f)  Solve  each  unknown  5  in  terms  of  5 2*  therefore 

5  4  =  5.4567833^ 

and 

=  13. 703306^2 

Step  4)  Skip  to  Rule  4  since  the  values  of 


and  are  unknown. 


T-3357 


a)  Write  the  equations  for  i.  and  *  . 

4  O 


6  =  X,n  °-4(5  +  j  ) 

4  l  ( '3  4* 


" 6  -  X2°  °‘6(45  +  V 


b)  Solve  each  equation  in  terms  of  X.. 


X1  " 


( £  +  s  ) 

\-3  -r  -4) 


(1/0-9162907) 


<S4> 


X2  = 


( *  5  +  *6> 


(V 


(1/0.5180256) 


c)  Replace  the  i's  with  known  values. 


X1  - 


(1  +  5 . 4567833 i^) 


( 5 . 4567833  5^ ) 


1.0913567 


X2  - 


(1  +  1 3 . 703306  $2 ) 


( 13.703306i2) 


1.9576152 
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d)  At  optimal ity  each  constraint  is  tight,  thus 
make  the  original  constraint  an  equality  and  solve 
it  in  terms  of  XjX2. 


X?X2  * 


.60 


e)  Substitute  the  equations  found  in  Step  4c  for 
XjX2  in  the  above  equation,  and  raise  each 
variable  to  the  appropriate  power. 


60 


(1  +  5. 4567833  ) 

(5.4567833S2) 


( 1 .0913567) (5) 


times 


(1  +  13. 7033065 2) 
(  13.703306?2) 


(1 .9576152) (7) 


Step  5)  The  above  equation  is  very  messy,  hence  all 
maximization  problems  regardless  of  the  number  of 
variables  will  be  solved  using  the  Brent-Dekker  method 
discussed  in  Chapter  3.  With  a  starting  bracket  (B,C) 
of  (0.001,1000)  the  value  for  ?2  Is  found  to  be 

5  2  =  0.260661 14 


Thus 


5  .  =  1 .422371  10 
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=  3.5719188 


Step  6)  Solve  for  Xj  and  X2  according  to  Rule  4  by 
replacing  the  S's  with  their  known  values. 


a)  Write  the  equations  for  5  .  and  S^. 

A  o 


‘4  ~  i  1  3  +  ’ 4 ' 


6  ~  *2  (  5  +  6 } 


b)  Substitute  in  the  known  values  for  the  S's. 


1.4223711  =  x“°-916290700( 1  +  1.4223711) 


3.5719188  =  x2°-05108256(.  +  3.5719188) 


c)  Solve  the  above  equations  for  X{  and  X2> 


Xj  =  60.913084 


X2  =  280.36967 


Step  7)  Use  the  above  values  to  solve  the  objective 


function. 


-1-  ~  f  *  S','  J*  sj'  m'J.  V «~  "Jv  .*L.  ."-V.  ■  .  s'-  r.'_  V.  s’-  "*-'**- 
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a)  Find  the  optimal  GP  values  of  and  by 
substituting  in  the  values  for  the  In  Xj. 


Nj  =  In  Xj 
N2  =  In  X2 


In  60.913084  =  4.109448 

In  280.36967  =  5.636109 


■j|A  b)  Therefore,  the  optimal  GP  value  of  the 

f.V/  objective  function  is 


maximum  reliability:  R  =  0.92195686 


Phase  3:  Integer  Solution 

Components  cannot  be  split  into  fractions,  therefore  an 
integer  solution  is  required.  This  phase  will  consist  of 
four  steps  that  will  ultimately  provide  the  maximum 
reliability  integer  component  solution. 

Step  1)  Round  each  Nj  down  to  the  nearest  integer,  and 
write  its  value. 


N 


1 


4 


N 


2 


5 


Step  2)  This  value  will  more  than  satisfy  the 
original  constraint;  however,  it  may  not  be  the  optimal 
combination  of  variables  that  both  satisfies  the 


'x  .  .  ■’’-V-V-  V  -"■f.y.  T-*.  tJlC ■> i"l  U 


constraint  and  has  the  maximum  reliability.  Determine 


the  cost  and  reliability. 

cost  =  $5(4)  +  $7(5)  =  $55.00  <  $60.00 


reliability  =  [l  -  0.4^)  (  1  -  0.6^]  =  0.8986307 

Step  3)  Determine  the  number  of  combinations  that 
must  be  compared  to  ensure  that  the  optimal  solution 
is  found. 

a)  Subtract  1  from  each  integer  value  found  in 
Step  1.  Determine  the  new  cost  for  each  element. 

cost  of  element  I  =  $5(3)  =  $15.00 

cost  of  element  2  =  $7(4)  =  $28.00 

b)  Use  the  fol lowing  formula  developed  by 
Steven  A.  Strauss  to  find  the  number  of 
combinations  to  be  checked  to  ensure  that  the 
maximum  reliability  is  found  subject  to  the 
maximum  al lowed  cost. 


+  combination  Nj  =  +  integer  of 


(budget  limitation  -  costnew) 


(cost  of  jth  element) 


-  i nteger  Nj  +  1 
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& 


where 


costnew  =  the  cost  of  each  element  found  in 
Step  3a,  plus  all  the  other  element  costs 
except  for  the  jth  element  (the  one  whose 
optimal  number  of  combinations  is  currently 
being  calculated) 

Therefore, 


+  combination  Nj  =  +  integer  of 


($60)  -  ($28) 

$5 


4  +  1  =  +  3 


+  combination  =  +  integer  of 


($60)  -  ($15) 

$7 


5  +  1  =  +  2 


c)  Find  all  combinations  for  each  variable  as 
determined  in  Step  3b.  The  combinations  are  also 
a  sensitivity  analysis  for  cost  and  reliability. 
Table  4.1  shows  the  combinations  of  Nj  and 


which  result  for  example  (1.2). 


Table  4.1  Maximum  Integer  Combinations 


N1 

N2 

COST 

(IN  DOLLARS) 

RELIABILITY 

2 

3 

31.00 

0.6585600 

2 

4 

38.00 

0.7311360 

2 

5 

45.00 

0.7746816 

2 

6 

52.00 

0.8008090 

2 

7 

59.00 

0.8164854 

3 

3 

36.00 

0.7338240 

3 

4 

43.00 

0.8146944 

3 

5 

50.00 

0.8632166 

3 

6 

57.00. 

0.8923300 

3 

7 

64.00 

0.9097980 

4 

3 

41.00 

0.7639296 

4 

4 

48.00 

0.8481178 

4 

5 

55.00. 

0.8986307 

4 

6 

62.00. 

0.9289384 

4 

7 

69.00 

0.9471230 

5 

3 

46.00 

0.7759718 

5 

4 

53.00 

0.8614871 

5 

5 

60 . 00. 

0.9127963 

5 

6 

67.00. 

0.9435818 

5 

7 

74.00 

0.9620531 

6 

3 

51.00 

0.7639296 

6 

4 

58.00 

0.8668348 

6 

5 

65.00. 

0.9184625 

6 

6 

72.00. 

0.9494391 

6 

7 

79.00 

0.9680251 

*  denotes  those 

combinations 

which  are  not 

feasible  since 

they  are  more 

than  the 

maximum  cost 

1  imitation.  The 

values  for  the 

combinations 

when  N .  is  I  and 

7  are  not 

included  as  the 

rel fabi 1  it ies 

are  too  low.  or 

the  cost  exceeds  the  budget  limitation. 
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Chapter  5 

SOLUTIONS  TO  APPLIED  SAMPLE 
RELIABILITY  PROBLEMS 

A  generalized  computer  algorithm,  as  described  in 
Chapters  3  and  4,  was  developed  for  a  particular  class  of 
nonlinear  reliabilty  design  problems  and  is  presented  in 
Appendix  B.  The  computer  algorithm  is  written  in  MICROSOFT 
Quickbasic  for  IBM-compatible  microcomputers. 

The  algorithm  is  built  around  5  areas:  (1) 
substitution  and  manipulation  of  terms  and  variables  to 
rewrite  the  original  problem  into  standard  GP  form;  (2)  the 
use  of  GP  to  start  solving  for  the  5  values;  (3)  the  use  of 
the  Brent-Dekker  method  to  solve  for  the  remaining  unkown  i 
values;  (4)  the  use  of  GP  to  solve  for  the  optimal 
non integer  answer;  and  (5)  solving  for  the  optimal  integer 
solution. 

The  type  of  reliability  problems  that  this  algorithm 
handles  (minimization  or  maximization)  will  always  remain 
at  one  DO.  no  matter  how  many  variables  are  added.  This  is 
due  to  the  substitution  and  manipulation  that  is  done. 
Whenever  another  variable  is  added,  one  more  constraint 
with  two  terms  (one  with  an  artificial  variable)  will  also 
be  added. 


,  *  u  r.  * \  *  . 


Currently  the  computer  algorithm  is  limited  to  four 
variables.  The  results  of  running  sample  problems  using 
the  computer  algorithm  are  addressed  below. 

Minimization  Problems 

Most  industries,  like  the  Department  of  Defense,  are 
limited  on  how  much  they  can  spend  to  either  purchase  new 
systems  or  upgrade  the  present  ones.  Electronic  equipment 
must  therefore  have  a  high  degree  of  reliability  with  a 
reasonable  cost  in  order  to  be  competitive  in  today's  world 
of  advanced  technology. 

The  following  problems  are  indicative  of  the 
reliability  constraints  imposed  on  civilian  contractors  for 
electronic  equipment  used  by  the  military.  The  first  two 
problems  were  developed  by  COL  Arbogast,  a  career  Signal 
Officer  currently  stationed  at  the  United  States  Military 
Academy.  The  third  problem  is  one  that  was  used  as  a 
sample  problem  in  a  published  reliability  text  book 
(Dhillon,  Balbir  1983).  Three  minimization  problems  are 
addressed  below  as  test  cases.  There  is  a  separate  test 
problem  depicting  each  number  of  variables  solvable  by  the 
computer  algorithm.  For  each  test  problem,  the  original 
problem  will  be  stated,  along  with  both  the  GP  and  integer 
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Two- Variable  Problem 

The  first  problem  depicts  a  system  that  is 
comprised  of  two  elements.  Table  5.1  provides  both 
the  optimal  GP  and  Integer  solutions.  The  original 
problem  is  shown  below. 


minimize  cost: 


subject  to: 


C  =  $40N  j  +  $60N2 


1' 


>  .98 


Table  5.1  Minimum  Two-Variable  Solutions 


Solution 

GP 

Integer 

Reliability 

0.979998 

0.989604 

Minimum  Cost 

$1 15.51 

$149.00 

N  j  Value 

1.228815 

2 

Va 1 ue 

0.9617115 

1 

Three-Variabl e  Problem 

This  problem  depicts  a  system  that  is  comprised 
of  three  elements.  Table  5.2  depicts  both  the  optimal 
GP  and  integer  solutions  for  the  problem.  The 
original  problem  Is  shown  below. 


fcA  > 
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m 


inimlze  cost:  C  =  $40Nj  +  $52N2  +  $60N^ 


subject  to: 


N  N 

1  -  0.01  1  -  0.006 


N. 


1  -  0.003 


>  0.97 


Table  5.2  Minimum  Three-Varlabl e  Solutions 


Solution 

GP 

Integer 

Rel iabi 1 Ity 

0.969840 

0.990919 

Minimum  Cost 

$133.84 

$192.00 

N j  Value 

1 . 020958 

2 

N2  Value 

0.888597 

1 

N3  Value 

0.779838 

1 

Foui — Variable  Problem 

This  problem  taken  from  Rel lab! 1 Ity  Engineering 
In  Systems  Design  and  Operation  ( pages  114-117),  by 
Balblr  Ohl I  Ion,  depicts  a  system  that  Is  comprised  of 
four  elements.  Table  5.3  depicts  both  the  optimal 
GP  and  integer  solutions  for  the  problem.  The  answers 
derived  by  the  computer  program  differ  slightly  from 
the  ones  given  In  the  test  book.  The  answers  that  are 


T-3357 


49 


different  in  the  text  book  are  in  parenthesis  by  the 
computer  answers.  It  appears  that  round  off  error  may 
cause  this  difference.  The  original  problem  is  shown 
below. 


minimize  cost:  C  =  $5Nj  +  SION^  +  $8N^  +  $2N4 


subject  to: 


N  N_| 

1  -  0.04  1  -  0.03 


1 

N„j 

1 

O 

• 

o 

M 

ui 

A 

i  -  .05 

>  0.98 


Table  5.3  Minimum  Four-Variabl e  Solutions 


Solution 

GP 

I nteger 

Rel iabi 1 Ity 

0.980302 

0.984330 

Minimum  Cost 

$95.54 

$100.00  ($98) 

N  j  Value 

5.61 1212 

6 

N2  Value 

3.924052 

4 

Va 1 ue 

3.859593 

4 

Va 1 ue 

8.332921 

9  (8) 

Max i m 1 zat i on  Prob 1 ems 

Cost  budgets  are  an  important  part  of  the  total  system 
planning,  especially  when  dealing  with  electronic  pieces  of 


r"wn 
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equipment.  Two  of  the  following  problems  were  developed  by 
COL  Arbogast,  while  the  third  was  taken  from  the  article 
published  by  Federowicz  and  Mumzudar.  Three  maximization 
problems  will  be  addressed  as  test  cases;  one  for  each 
number  of  variables  solvable  by  the  computer  algorithm. 

Two- Variable  Problem 

The  first  problem  is  one  depicting  a  system  that 
is  comprised  of  two  elements.  Table  5.4  provides  both 
the  optimal  GP  and  integer  solutions.  The  original 
problem  is  shown  below. 


N 

maximize  reliability:  R  =  | 1-  0.02 
subject  to:  $40Nj  +  $60N2  <  $200 


N. 


1  -  0.01 


Table  5.4  Maximum  Two-Variable  Solutions 


Solution 


GP  Integer 


Reliability 
Minimum  Cost 
N j  Value 
N_  Value 


0.999635 
$200.00 
2.233251 
1 .844499 


0.999500 

$200.00 

2 

2 


This  problem  depicts  a  system  that  is  comprised 


of  three  elements.  Table  5.5  depicts  both  the  optimal 
GP  and  integer  solutions  for  the  problem.  The  original 
problem  is  shown  below. 


(  Nl] 

1  N2 

max i m i ze  re  1 i ab i 1 i t  v :  R  =  1 1  -  0.01 

i 

>1  -  0.004  *■] 

Il  -  0.003  ] 

subject  to:  $40Nj  +  $52N^  +  $60N^  <  $250 

Table  5.5  Maximum  Three- Var iable  Solutions 


So  I ut i on 


GP  Integer 


Reliability 
Minimum  Cost 
N  Value 
Va 1 ue 
Value 


0.999546 
$250.00 
1.928819 
1.59408 
1 . 49925 


0.9968843 

$244.00 

2 

2 

1 
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Four-Var i ab 1 e  Problem 

This  problem  taken  from  the  art i c I e  published  bv 
Federowi cz  and  Mazumdar,  depicts  a  svstem  that  is 
comprised  of  four  elements.  Table  5.3  depicts  both  the 
optimal  GP  and  the  integer  solutions  for  the  problem. 
The  answers  derived  by  the  comput  •*  program  differ 
sliahtlv  from  the  ones  given  in  the  test  article. 
Roundoff  error  and  the  way  combinations  are  checked 
appear  to  be  the  reason  for  the  computer  not  attaining 
the  correct  optimal  answer.  The  answers  that  differ 
are  noted  with  the  article  answer  in  parentheses.  The 
original  problem  is  shown  below. 


maximize  reliability:  R  = 


N, 

N_ 

1  -  0.20 

2 

!  -  0.30 

N. 


1  -  0.25 


N. 


1  -  0.15 


subject  to: 


$  1 . 20N |  +  $2.3N2 


+  $3 . 40N  +  $4 . 50N  <  *47 

3  4  — 


T-3357 


Table  5.6  Maximum  Foui — Variable  Solutions 


Solution 

GP 

1 nteger 

Reliability 

0.993107 

0.988735  (0.99) 

Minimum  Cost 

$47.00 

$43.40  ($45.80) 

Nt  Value 

4.522890 

4  (6) 

Value 

5.265519 

5 

Value 

4.393139 

4 

N4  Value 

3.227793 

3 
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Chapter  6 
CONCLUSIONS 

Technology  Is  an  ever  changing  area  of  science.  As 
equipment  becomes  more  sophisticated,  the  need  for  reliable 
and  cost  effective  elements  Increases.  A  reliability  model 
that  is  quick,  accurate,  and  efficient  is  an  important  part 
of  resource  management.  It  is  evident  that  the  generalized 
computer  algorithm  can  be  a  useful  microeconomic  tool  in 
terms  of  convenience,  speed  of  results,  and  reduction  of 
computational  error  in  calculating  the  optimal  combinations 
of  system  elements. 

The  strength  of  this  computer  algorithm  is  that  it 
provides  a  method  for  solving  difficult  problems  that  are 
currently  either  guessed  at,  or  require  an  excessive  amount 
of  time  to  solve  either  by  hand  or  computer.  The  program 
is  easy  to  use  and  solves  each  problem  very  quickly.  The 
algorithm  has  a  slight  weakness.  When  the  costs  for  each 
element  are  relatively  small,  with  little  difference,  and 
the  failure  rates  are  also  fairly  close,  the  computer 
sometimes  does  not  recognize  a  difference  between  adding  or 
subtracting  one  element.  This  appears  to  happen 
occassional ly  when  there  are  four  variables. 


-■  <*  . , . , 
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One  area  not  ful ly  explored  was  that  of  sensitivity 
analysis.  The  computer  algorithm  presently  checks 
different  combinations  of  the  variables,  but  only  with 
their  present  cost  and  failure  rate  values.  The  algorithm 
should  be  expanded  to  include  a  method  that  would  determine 
how  sensitive  each  element  is  to  its  own  cost  and  failure, 
as  wel 1  as  that  of  the  total  system. 

Another  important  area  for  further  research  is  that  of 
expanding  the  computer  algorithm  to  include  more  than  four 
variables.  This  means  developing  a  method  to  determine  the 
appropriate  function  to  be  solved  when  searching  for 
The  subrout i ne  that  searches  for  the  unknown  5  va 1 ues 
should  also  be  refined.  The  Brent-Oekker  method  is  very 
sensitive  when  converging  on  the  positive  root.  When  the 
method  finds  a  positive  root,  the  convergence  factor 
becomes  less  concise  as  the  number  of  variables  increases. 
It  is  also  noted  that  as  the  required  minimum  reliabilities 
become  more  precise,  the  less  concise  the  criteria  for  the 
convergence  factor  can  be .  A  1 ower  convergence  factor  i s 
therefore  necessary,  otherwise,  the  computer  tends  to 
osci 1 1 iate  around  the  positive  root. 

The  program  could  also  be  modified  to  solve  problems 
with  additional  constraints.  For  Instance,  there  may  be  a 
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size  capacity  which  will  require  a  constraint  limiting  the 
total  number  of  components. 

Another  area  for  further  research  is  that  of  improving 
the  manipulation  of  the  original  problem  to  determine 
whether  a  different  method  will  reduce  the  problem  to  a 
zero  DD,  GP  problem.  If  the  problem  could  be  reduced  to 
zero  DD  it  would  cut  the  solving  time  in  half  as  the  Brent- 
Dekker  method  could  then  be  eliminated  from  the  program. 

The  concept  of  using  manipulation  and  substitution  to 
form  messy  problems  into  ones  solvable  by  GP,  and  then 
using  a  nonlinear  method  such  as  Brent-Dekker  to  solve  the 
Inherent  one  DD  mathematical  problem  significantly  broadens 
the  range  of  possible  applications  of  GP.  This  in  itself 
is  an  important  area  for  further  research. 

The  main  contribution  of  this  thesis  is  that  it  will 
provide  a  simple-to-use,  generalized  computer  algorithm  for 
solving  a  large  class  of  nonlinear  equations.  In  addition, 
the  algorithm  guarantees  that  the  GP  solution  will  be 
globally  optimal.  The  integer  solution  cannot  guarantee 
global  optimality,  however,  in  most  cases  it  will  provide 
the  optimum  local  integer  combination. 
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Appendix  A 


FOUR  RULES  FOR  GEOMETRIC  PROGRAMMING 

Rule  1 :  The  form  of  the  optimal  solution  of  any  posynomtal 

GP  prob 1 em  is: 

Value  of  the  Objective  Function  = 

A 

(coef.  of  first  term  In  obj.  f unct 1 on/ 5  ) 1  X 

...  X 

J  1  ast 

(coef.  of  Fast  term  fn  obj.  funct fon/5 ,ast >  X 

(coef.  of  first  term  in  constralnt/5 )ast+1 > ,ast+1  x 
...  X 

(coef.  of  last  term  In  constraint/ !  last+m  x 

I ast+m 

(sum  of  S's  In  constraint) 
(sum  of  - ' s  in  constraint) 


Rule  2:  The  exponent  matrix  Is  constructed  in  the  following 
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Rule  2A:  The  sum  of  contributions  to  cost  in  the 
objective  function  =  1. 


51  +  S2  + 


+  5 


last 


1 


Rule  2B:  For  each  primal  variable  the  equations  in 
the  exponent  matrix  are: 


(power 

of 

var i abl e 

j 

in  term 

1  ) 

X  S,  ♦ 

(power 

of 

var i ab 1 e 

j 

i n  term 

2) 

x  «2  - 

(power  of  variable  J  in  last  term)  X 


last+m 


=  0 


Rule  3 :  At  opt  I ma 1 i ty 

Value  of  the  Objective  Function  = 

(first  term  in  obj.  functlon/5j)  = 
(second  term  in  obj.  function/ S^)  = 


(last  term  in  obj.  funct ion/ 5 (  ^ > 

Rule  4 ;  At  optimality,  for  each  constraint 

i.  =  ( Jth  term  in  constraint)  X 


(sum  of  i's  in  constraint) 
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Append i x  B 

COMPUTER  ALGORITHM  FOR  SOLVING 
RELIABILITY  PROBLEMS  WITH  GP 
USING  MICROSOFT  QUICKBASIC 

The  reliability  program  developed  for  this  thesis  is 
designed  to  run  on  an  I BM-compat ibl e  personal  computer 
using  a  software  package  called  Microsoft  Ouickbasic.  To 
use  the  Microsoft  Ouickbasic  package,  the  user  must  have 
access  to  a  computer  with  at  least  384k  memory.  The 
computer  program  could  easily  be  modified  and  written  in 
the  Basic  language  of  the  user's  personal  computer  if 
M 1 crosof t  Qu 1 ckbas 1 c  is  not  ava liable. 

To  run  the  program,  the  program  Is  first  read  into  the 
computer  and  an  executable  file  is  made.  The  program  can 
be  then  be  executed  by  typing  the  file  name  (i.e.,  RELY)  at 
the  computer  A  prompt. 

The  program  is  designed  to  be  user  interactive.  Once 
the  program  is  loaded,  the  computer  will  begin  by  asking 
the  user  for  input.  Appendix  C  depicts  a  sample  computer 
run. 

The  reliability  computer  algorithm  using  Microsoft 
Ouickbasic  is  listed  below. 


wJv.iit 


»*.v 

-w 

-w 

>.v 

Ui.&> 


SR 

❖ 


T-3357 


WELCOME  TO  CPT  OATNEY'S  RELIABILITY  SOLVING 
PROGRAM.  THIS  PROGRAM  WILL  EITHER  MINIMIZE 
COST  SUBJECT  TO  A  REQUIRED  RELIABILITY  OR 
MAXIMIZE  RELIABILITY  SUBJECT  TO  A  BUDGET 
LIMITATION  USING  GP.  THE  PROGRAM  WILL  HANDLE 
FROM  2  TO  4  VARIABLES.  THE  FINAL  ANSWER,  IN 
INTEGER  FORM,  WILL  INCLUDE:  THE  REQUIRED 
QUANTITY  FOR  EACH  VARIABLE  AND  EITHER  THE 
MINIMUM  COST  OR  MAXIMUM  RELIABILITY  OBTAINED 
FOR  THE  DEFINED  PROBLEM. 


OPTION  BASE  1 
DEFINT  I-L.N 
DEFDBL  A-D.F.M.P 

DIM  COF (10) ,PROB( tO).XX(IO) ,YY( 10) .DELTA (20) 
DIM  NO( 10) ,PROBNEW( 10) ,DIFFER( 10) . JCO( 10) 

DIM  POWER (10), COSTNEW (10) 


CLS 

GOSUB  START I NP 


IF  YN$="N"  THEN  GOSUB  MAXIM  ELSE  GOSUB  MINIM 


GOSUB  ROUNDOFF 
GOSUB  VERIFY 
GOSUB  PRINTOUT 
GOTO  9999 


START I NP: 


THIS  ROUTINE  INITIALIZES  THE  PROGRAM  **** 

.  BY  ASKING  THE  USER  TO  INPUT  ******* 

. *  THE  ORIGINAL  PROBLEM  DATA  . 


10  INPUT;  "IS  THIS  PROBLEM  A  MINIMIZATION?  (Y/N)  ",Y 

IF  YN$< >"Y"  AND  YN$<>"N"  THEN  GOTO  10 
PRINT  "" 

15  INPUT;  "HOW  MANY  VARIABLES  ARE  THERE?  (2-4)  ",NVA 
PRINT  "" 

16  PRINT  "YOU  HAVE  INPUT  " , NVAR 
PRINT  "AS  THE  NUMBER  OF  VARIABLES" 

INPUT;  "IS  THIS  VALUE  ACCEPTABLE?  (Y/N)  " , YYNN$ 
PRINT  "" 

IF  YYNN$="N"  THEN  GOTO  15 
IF  YYNN$< >"Y"  THEN  GOTO  16 

FOR  J=1  TO  NVAR 

20  CLS 

LOCATE  5, 12 

PRINT  "YOU  ARE  INPUT I NG  VARIABLE  #".J 

INPUT;  "WHAT  IS  THE  COST  COEFFICIENT?  ",COF(J) 

PRINT  "" 

INPUT;  "WHAT  IS  THE  ASSOCIATED  PROBABILITY  OF 
FAILURE?  " * PROB ( J ) 

PRINT  "" 

PRINT  "YOU  HAVE  INPUT  ";COF(J);"  AND  ";PROB(J) 
PRINT  "AS  THE  COST  AND  FAILURE  PROBABILITY" 

30  LOCATE  15,12 

INPUT;  "ARE  THESE  VALUES  ACCEPTABLE?  (Y/N)  ",Y 
PRINT  "" 

IF  YNN$="N"  THEN  GOTO  20 
IF  YNN$<>"Y"  THEN  GOTO  30 
NEXT  J 
CLS 

40  IF  YN$="Y"  THEN 

INPUT;  "WHAT  IS  THE  REQUIRED  RELIABILITY?  ".RE 
PRINT  "" 

PRINT  "YOU  HAVE  INPUT  " ; RELY 
PRINT  "AS  THE  REQUIRED  RELIABILITY" 


T-3357 


63 


50  LOCATE  15,12 

INPUT;  "IS  THIS  VALUE  ACCEPTABLE?  (Y/N)  ' , YNN$ 
PRINT  "" 

IF  YNN$="N"  THEN  GOTO  40 
IF  YNN$< >"Y"  THEN  GOTO  50 


60  INPUT;  "WHAT  IS  THE  BUDGET  LIMITATION?  " , BUDG 

PRINT  "" 

PRINT  "YOU  HAVE  INPUT  " ; BUDG 
PRINT  "AS  THE  BUDGET  LIMITATION" 

70  LOCATE  15,12 

INPUT;  "IS  THIS  VALUE  ACCEPTABLE?  (Y/N)  " , YNN$ 
IF  YNN$="N"  THEN  GOTO  60 
IF  YNN$<  >"Y"  THEN  GOTO  70 
END  IF 
RETURN 


MAXIM: 


THIS  SUBROUTINE  PROVIDES  THE 
*  EQUATIONS  NEEDED  TO  SOLVE 
*****  FOR  THE  5  VALUES  IN  ** 
****  MAXIMIZATION  PROBLEMS  * 


DELTA ( 1 )=1 

FOR  J  =  1  TO  NVAR 

COFF ( J ) =COF ( J ) /BUDG 
POWER(J)  =  -COFF ( J ) /LOG ( PROB ( J ) ) 
DELTA (2*J+1 )=1 
NEXT  J 
GOSUB  DFUN 
GOSUB  SECSOLVE 
GOSUB  XMAXSOLVE 
RETURN 


**  THIS  SUBROUTINE  PROVIDES  THE  . 

****  EQUATIONS  NEEDED  TO  SOLVE  ****** 
********  FOR  THE  *  VALUES  IN  ******** 
*******  MINIMIZATION  PROBLEMS  ******* 
****  IT  ALSO  SOLVES  TWO-VARIABLE  **** 
PROBLEMS  USING  THE  QUADRATIC  EQUATION 


CON  =  -1 
DELTA! 1)  =  1 
AA  =  ( I -RELY) 

BB  =  0  :  CC  =  0  :  DD  =  0 
FOR  J  =  1  TO  NVAR 

DELTA(2*J+1)  =  -COF { J ) /LOG ( PROB ( J ) ) 

BB  =  BB-DELTA(2*J+1>*RELY 
CON  =  CON*DELTA( 2*J+1 ) *RELY 
FOR  IJ  =  J+l  TO  NVAR 

IF  I J>NVAR  THEN  GOTO  100 

CC  »  CC-DELTA(2*J+1 ) *DELTA( 2* I J+l )*RELY 

FOR  IJK  =  IJ+1  TO  NVAR 

IF  IJtONVAR  THEN  GOTO  100 
00  *  00-DELTA (2*J+1) -DELTA (2* I J+l )* 
DELTA ( 2* I JK+ 1 ) *RELY 

NEXT  IJK 
NEXT  IJ 

NEXT  J 

IF  NVAR=2  THEN 

02  »  (-BB+SQR(BB“2-4*AA*CON) )/(2*AA) 

GOSUB  XMINSOLVE 
RETURN 
END  IF 

IF  NVAR>2  THEN  GOSUB  DFUN 
GOSUB  SECSOLVE 
GOSUB  XMINSOLVE 
RETURN 
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DFUN : 


THIS  SUBROUTINE  DEFINES  THE  FUNCTION 
*  TO  BE  USED  WHEN  SOLVING  FOR  THE  *** 

****  REAL  POSITIVE  ROOT  OF  THE  . * 

****  UNKNOWN  5  VALUES  IN  EITHER  . 

MAXIMIZATION  OR  MINIMIZATION  PROBLEMS 


DEF  FNFUNC  (D2) 

IF  YN$="Y"  THEN 

FNFUNC  =  AA*D2"'NVAR+BB,D2'>  ( NVAR-  1  )+CC*D2~  ( NVAR-2 ) 
+DD*D2~ ( NVAR- 3 ) * ( NVAR- 3 ) +CON 

ELSE 

DUM=  1 

FOR  J  =  1  TO  NVAR 

DUM  =  DUM* ( ( ( 1 +POWER ( J ) *D2 ) / 

(POWER!  J)*D2)  KPOWER(J)  ) 

NEXT  J 

FNFUNC  =  -EXP( 1 )+OUM 
END  IF 
END  DEF 
RETURN 


SECSOLVE : 


THIS  SUBROUTINE  FINDS  THE  POSITIVE 
***  REAL  ROOT  FOR  THE  UNKNOWN  5 'S  * 
THE  SUBROUTINE  BEGINS  BY  REQUESTING 
•***  AN  INITIAL  STARTING  BRACKET  ** 


220  PRINT  "PLEASE  PROVIDE  A  STARTING  BRACKET  [B,C] 
SEPARATED  BY  A  COMMA  " 

INPUT  "WHICH  YOU  BELIEVE  CONTAINS  A  ROOT  ";B,C 
IF  B<0  OR  C<0  OR  (B=0  AND  C=0)  THEN  GOTO  220 


240  PRINT  "THANK  YOU,  YOUR  STARTING  BRACKET  IS  =  ",B,C 
D=C : GOSUB  900 
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. . *****  MAiN  PROGRAM  STARTS  HERE  ******** 

'***  QUICK  CHECK:  IS  FA  AND  FB  OF  OPPOSITE  SIGNS? 

A=C 

GOSUB  760 
GOSUB  690 

IF  FA*FB  >  0  THEN  GOTO  1350 
CONVERG=lO~(-8+(2*(NVAR-2) )) 

IF  ABS(FB*FD)<CONVERG  THEN 
IF  ABS ( FD ) < ABS ( FB )  THEN 
B=D 
FB=FD 
END  IF 
END  IF 

IF  ABS ( FB ) <CONVERG  THEN  1200 
GOSUB  830 

IF  ABS(FC) =>ABS(FB)  THEN  GOTO  480 
A=8 

GOSUB  760 
B=C 

GOSUB  690 
C=A 

GOSUB  830 
GOTO  350 


PRINT  CURRENT  VALUES  OF  B,A,FB,FA  **< 


PRINT 

D=(FB*A-FA*B)/(FB-FA) 

PRINT 

GOSUB  900 

A=B 

PRINT 

GOSUB  760 

PRINT 

M=(B+C)/2 

GOSUB  920 

GOTO  970 


FB=FNFUNC(B) 

PRINT 

RETURN 


CALCULATE  FB 
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CALCULATE  FA 


FA=FNFUNC( A) 

PRINT 

RETURN 


CALCULATE  FC 


FC=FNFUNC(C) 

PRINT 

RETURN 


CALCULATE  FO 


FD=FNFUNC(D) 

PRINT 

RETURN 


CALCULATE  FM 


FM=FNFUNC(M) 

PRINT 

RETURN 


PRINT  CURRENT  VALUES  OF  D.B.M.AND  C 


PRINT 

PRINT  -THE  CURRENT  VALUES  OF  D,B,M,C  ARE  = 

FXX=FD*F8 

FYY=FD*FM 

FZZ=FC*FD 

IF  FYY<=0  THEN 

IF  ABS(FD) <ABS(FM)  THEN 
B=D :FB=FD 
C=M:FC=FH 
ELSE 

C=0:FC=FD 
B=M : FB=FM 
ENO  IF 
A=C :FA=FC 
GOTO  330 


,D,B,M,C 


iw.if 


ELSE 

IF  FZZ<=0  THEN 
B=C:C=A: A=B 
FB=FC:FC=FA:FA=FB 
GOTO  1050 
END  IF 
END  IF 
PRINT 

1050  IF  (D>B  AND  D<M)  AND  FXX<0  THEN  GOTO  1060 
IF  (D<B  AND  D>M)  AND  FXX<0  THEN  GOTO  1060 
GOTO  1 1 00 

1060  B=D 

FB=FD 
GOTO  1 1 20 

1100  B=M 

FB=FM 

1120  IF  (FB*FC)=>0  THEN  C=A:GOSUB  830 
GOTO  330 


'•»••••••••••••••••  print  RESULTS  **************** 

1200  CLS 

PRINT 

PRINT  "THE  OPTIMAL  VALUE  OF  THE  UNKNOWN  f  IS  =  ".B 

PRINT 

D2=B 

GOTO  1400 


'***********  TELL  USER  BRACKET  IS  NO  GOOD . . 

1350  PRINT 

PRINT  "ROOT  NOT  BRACKETED  OR  SECANT  METHOD  WILL  NOT 
FIND  ROOT,  TRY  AGAIN" 

GOTO  220 


USER  IS  FINISHED  WITH  SUBROUTINE 


1 400  RETURN 
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XMAXSOLVE 


THIS  SUBROUTINE  SOLVES  FOR  THE  X  VALUES  **** 
► .  IN  MAXIMIZATION  PROBLEMS  ************* 


FOR  J  =  1  TO  NVAR 

DELTA(2*J+2)=D2*POWER( J) 

XX(J)  =  ( LOG ( DELTA! 2* J+2 ) / (DELTA ( 2* J+ I ) 
+OELTA ( 2* J+2 ) ) ) )/LOG(PROB(J) ) 

NEXT  J 
RETURN 


XM INSOL VE : 


THIS  SUBROUTINE  SOLVES  FOR  THE  X  VALUES 
»••••»  jN  MINIMIZATION  PROBLEMS  . 


FOR  J  =  I  TO  NVAR 
DELTA! 2* J ) =D2 

XX(J)  =  ( LOG (DELTA! 2* J+ 1 ) / ( DELTA ( 2* J ) 
+DELTA ( 2* J+ 1 ) ) ) ) /LOG ( PROB ( J ) ) 

Ni_XT  J 
RETURN 


ROUNDOFF : 


THIS  SUBROUTINE  ROUNDS  THE  GP 
*•  SOLUTION  TO  INTEGER  VALUES  * 


FOR  J  =  I  TO  NVAR 
YY ( J )  =  I  NT ( XX ( J ) } 

IF  INT(XX(J) )<XX( J)  AND  YN$=' 
YY ( J )  =  I  NT ( XX ( J ) )  +  l 

NEXT  J 
RETURN 


THEN 


VERIFY: 


THIS  SUBROUTINE  FINOS  AND  VERIFIES 
**  THE  OPTIMAL  INTEGER  SOLUTION  * 


PROBOPT  = 1 :COSTOPT=0 
FOR  J  =  1  TO  NVAR 

COSTOPT=COSTOPT+COF ( J ) *  Y  Y ( J ) 

PROBOPT = PROBOPT* ( 1 -PROB ( J ) ~YY ( J ) ) 

NEXT  J 

IF  YN$="N”  THEN 

FOR  J  =  1  TO  NVAR 
NO( J ) =YY ( J ) -  1 
NEXT  J 

FOR  J  =  1  TO  NVAR 
COSTNEW ( J ) =0 
FOR  K  =  1  TO  NVAR 
IF  JOK  THEN 

COSTNEW ( J ) =COSTNE W ( J  >  +  ( COF ( J ) *  NO ( K ) ) 

NEXT  K 

DIFFER! J)  =  I  NT ( ( BUDG-COSTNEW ( J ) ) / 

COF{ J) ) -YY ( J ) + 1 

NEXT  J 
END  IF 

IF  YN$="Y”  THEN 

FOR  J  =  1  TO  NVAR 
NO( J )=YY(J)+I 
NEXT  J 

FOR  J  =  1  TO  NVAR 
PROBNEW! J}=1 
FOR  K  =  1  TO  NVAR 

IF  JOK  THEN  PROBNEW  (  J )  =PROBNEW  ( J )  * 

( l-PROB(K) 'NO(K) ) 

NEXT  K 

OIFFER(J) =YY ( J ) -  I NT ( (LOG( 1 - ( RELY/ 

PROBNEW! J) ) ) ) / ! LOG ( PROB ! J ) ) ) ) 

NEXT  J 
END  IF 


* 

*» 

*1 
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TEST  FOR  CONSTRAINT  VIOLATION  AND 
**  A  BETTER  OBJECTIVE  FUNCTION  ** 


FOR  K  =  1  TO  NVAR:  NO(K)=YY(K):  NEXT  K 


DETERMINE  THE  NUMBER  OF  COMBINATIONS 
*****  TO  BE  EXAMINED  TO  FIND  ******* 
******  THE  OPTIMAL  SOLUTION  ******** 


FOR  J1  =  YY ( 1 ) -D I FFER ( I )  TO  YY ( 1 )+DIFFER( I ) 

JCO( 1 )=J1 

IF  JKI  THEN  GOTO  9900 

FOR  J2  =  YY(2)-DIFFER(2)  TO  YY(2)+DIFFER(2) 

JCO( 2 ) =J2 

IF  J2< 1  THEN  GOTO  9800 

IF  NVAR=2  THEN  GOSUB  COMPARE:  GOTO  9800 
FOR  J3  =  YY ( 3 ) -D I FFER ( 3 )  TO  YY(3)+DIFFER(3) 
JCO(3)=J3 

IF  J3< 1  THEN  GOTO  9700 

IF  NVAR=3  THEN  GOSUB  COMPARE:  GOTO  9700 
FOR  J4=YY(4)-DIFFER(4)  TO  YY(4)+DIFFER(4) 
JCO(4)=J4 

IF  J4< 1  THEN  GOTO  9600 
GOSUB  COMPARE 
9600  NEXT  J4 

9700  NEXT  J3 

9800  NEXT  J2 

9900  NEXT  J 1 
RETURN 


COMPARE : 


THIS  SUBROUTINE  FINDS  THE 
**  OPTIMAL  COMBINATION  ** 


IF  YN$="N"  THEN 
COSTCALC=0 
FOR  JJ  =  1  TO  NVAR 

COSTCALC=COSTCALC+COF ( J J ) * JCO ( J J ) 
NEXT  JJ 


T-3357 


IF  COSTCALOBUDG  THEN  RETURN 

PROBCALC= 1 

FOR  JJ  =  1  TO  NVAR 


PROBCALC=PROBCALC* ( l-PROB(JJ) ~JCO< JJ) ) 
NEXT  JJ 

IF  PROBC ALC < PROBOPT  THEN  RETURN 
PROBOPT  =  PROBCALC 
FOR  JJ  =  I  TO  NVAR 
NO( J J ) = JCO ( J J ) 

NEXT  JJ 
ELSE 

PROBCALC= 1 

FOR  JJ  =  1  TO  NVAR 

PROBCALC=PROBCALC* ( 1 -PROB ( J J ) * JCO< JJ ) ) 
NEXT  JJ 

IF  PROBCALC < RELY  THEN  RETURN 

COSTCALC=0 

FOR  JJ  *  1  TO  NVAR 

COSTCALC=COSTCALC+COF ( J J ) « JCO ( J J ) 

NEXT  JJ 

IF  COST C AL C < COSTOPT  THEN 
COSTOPT  =  COSTCALC 
FOR  JJ  *  I  TO  NVAR 
NO( JJ)=JCO( JJ) 

NEXT  JJ 
END  IF 
ENO  IF 
RETURN 


PRINTOUT: 


THIS  SUBROUTINE  TELLS  THE  USER  " 
WHAT  THE  OPTIMAL  SOLUTION  IS  *" 


PRINT  "THE  ORIGINAL  PROBLEM  WAS:" 

IF  YN$="N"  THEN 

PRINT  "MAXIMIZE  RELIABILITY": 

FOR  JJ  =  I  TO  NVAR- I 

PRINT  " { 1 - ( " ; PROB( JJ );"~N(";JJ "))*"; 
NEXT  JJ 

PRINT  "( 1 - ( " ; PROB ( NVAR ) ; "“N( " ; NVAR" ) ) " ; 
PRINT  "  S.T.  "; 

FOR  JJ  =  1  TO  NVAR- I 

PRINT  COF(JJ) ;"*N(";JJ:") 
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39 

m 


•art 


a 


•A 

c 


9997 


NEXT  JJ 

PRINT  COF(NVAR) ; "*N( " ;NVAR; " )  <=  ";BUOG 
PRINT 
ELSE 

PRINT  "MINIMIZE  COST"; 

FOR  JJ  =  1  TO  NVAR-1 

PRINT  COF (JJ) ;"N("; JJ;")  +"; 

NEXT  JJ 

PRINT  COF(NVAR) ; "*N( " ; NVAR; " ) " 

PRINT  "  S.T. 

FOR  JJ  =  I  TO  NVAR- I 

PRINT  " ( 1 - ( " ; PROB (JJ) ;"~N(";JJ") )*"; 

NEXT  JJ 

PRINT  "( l-(";PROB(NVAR) ; "~N ( " ; NVAR" ) )  =>  " ; RELY 
PRINT 
ENOIF 

COSTOLD=0  ;  COSTCALOPT=0  :  PROBOLD=l  :  PROBCALOPT= 1 
PRINT  "THE  GP  SOLUTION  FOR  THE  PROBLEM  HAS  VALUES 
OF" 

FOR  J  =  1  TO  NVAR 

PRINT  "  N(";J;")=";XX(J) 

COSTOLD=COSTOLO+XX ( J ) *COF ( J ) 
COSTCALOPT=COSTCALOPT+NO ( J ) *COF ( J ) 
PROBOLO=PROBOLD* ( 1 -PROB ( J ) “XX ( J ) ) 
PROBCALOPT=PROBC ALOPT * ( I - PROB ( J ) “ NO ( J ) ) 

NEXT  J 

PRINT  "  WITH  AN  OPTIMAL  COST  OF  "; COS TOLD 
PRINT  "  AND  A  RELIABILITY  OF  "jPROBOLD 
PRINT 

PRINT  "THE  OPTIMAL  INTEGER  SOLUTION  HAS  VALUES  OF  " 
FOR  J  =  1  TO  NVAR 

PRINT  "  N(";J;")=" ; NO( J ) 

NEXT  J 

PRINT  "  WITH  AN  OPTIMAL  COST  OF  ";COSTCALOPT 
PRINT  "  AND  A  RELIABILITY  OF  ";PROBCALOPT 

INPUT  "DO  YOU  WANT  A  HARD  COPY  OF  THIS  RESULT? 

(Y/N)  " , NY$ 

IF  NY$="N"  THEN  GOTO  9998 
IF  NY$< >"Y"  THEN  GOTO  9997 

LPRINT 

LPRINT  "THE  OPTIMAL  VALUE  OF  IS  =  ",B 
LPRINT 

LPRINT  "THE  ORIGINAL  PROBLEM  WAS:" 

IF  YN$= "N"  THEN 

LPRINT  "MAXIMIZE  RELIABILITY"; 

FOR  JJ  =  I  TO  NVAR- I 


•  ww 


r? i 


PRINT  "( l-(";PROB(JJ) ;"~N(";JJ") )*"; 

NEXT  JJ 

LPRINT  "( l-<";PROB(NVAR) ; "~N ( " ; NVAR" ) )"; 

LPRINT  "  S.T.  "; 

FOR  JJ  =  1  TO  NVAR-1 

LPRINT  COF(JJ) ;"*N("; JJ;")  +"; 

NEXT  JJ 

LPRINT  COF ( N VAR ) ; " * N ( " ; N VAR ; " )  <=  " ; BUOG 
LPRINT 
ELSE 

LPRINT  "MINIMIZE  COST"; 

FOR  JJ  =  1  TO  NVAR-1 

LPRINT  COF(JJ) ;"*N("; JJ;")  +"; 

NEXT  JJ 

LPRINT  COF ( NVAR ) ;"*N("; N VAR ; " ) " 

LPRINT  "  S.T.  "; 

FOR  JJ  =  1  TO  NVAR- 1 

LPRINT  " ( 1 - ( " ;PROB( J J );"~N(";JJ "))*"; 

NEXT  JJ 

LPRINT  "( l-(";PROB(NVAR) ; "~N( " ;NVAR" ) )  =>  " ; RELY 
LPRINT 
END  IF 

LPRINT  "THE  GP  SOLUTION  FOR  THE  PROBLEM  HAS  VALUES 
OF" 

FOR  J  *  l  TO  NVAR 

LPRINT  "  N( " ; J; ")=";XX( J) 

NEXT  J 

LPRINT  "  WITH  AN  OPTIMAL  COST  OF  ";COSTOLD 
LPRINT  "  AND  A  RELIABILITY  OF  "; PROBOLD 
LPRINT 

LPRINT  "THE  OPTIMAL  INTEGER  SOLUTION  HAS  VALUES 
OF  " 

FOR  J  =  I  TO  NVAR 

LPRINT  "  N ( " ; J ; " ) =" ; NO( J ) 

NEXT  J 

LPRINT  "  WITH  AN  OPTIMAL  COST  OF  ";COSTCALOPT 
LPRINT  "  AND  A  RELIABILITY  OF  ";PROBCALOPT 

9998  RETURN 

9999  INPUT;  "DO  YOU  WISH  TO  DO  ANOTHER  PROBLEM?  (Y/N)  ", 

NY$ 

IF  NY$="Y"  THEN  GOTO  I 
IF  NY$< >"N"  THEN  GOTO  9999 
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Appendix  C 


SAMPLE  COMPUTER  RUN 


njter  Questions 


Input 


IS  THIS  PROBLEM  A  MINIMIZATION?  (Y/N) 

HOW  MANY  VARIABLES  ARE  THERE?  (2-4) 

YOU  HAVE  INPUT  3  AS  THE  NUMBER  OF  VARIABLES 
IS  THIS  VALUE  ACCEPTABLE?  (Y/N) 


YOU  ARE  INPUTTING  VARIABLE  1 
WHAT  IS  THE  COST  COEFFICIENT? 

WHAT  IS  THE  ASSSOCIATED  PROBABILITY  OF  FAILURE? 
YOU  HAVE  INPUT  5  AND  .4 
AS  THE  COST  AND  FAILURE  PROBABILITY 
ARE  THESE  VALUES  ACCEPTABLE?  (Y/N) 
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YOU  ARE  INPUTTING  VARIABLE  2 
WHAT  IS  THE  COST  COEFFICIENT? 

WHAT  IS  THE  ASSSOCIATED  PROBABILITY  OF  FAILURE? 
YOU  HAVE  INPUT  7  AND  .6 
AS  THE  COST  AND  FAILURE  PROBABILITY 
ARE  THESE  VALUES  ACCEPTABLE?  (Y/N) 


YOU  ARE  INPUTTING  VARIABLE  3 
WHAT  IS  THE  COST  COEFFICIENT? 

WHAT  IS  THE  ASSSOCIATED  PROBABILITY  OF  FAILURE? 
YOU  HAVE  INPUT  6  AND  .5 
AS  THE  COST  AND  FAILURE  PROBABILITY 
ARE  THESE  VALUES  ACCEPTABLE?  (Y/N) 


WHAT  IS  THE  REQUIRED  RELIABILITY? 
YOU  HAVE  INPUT  .90 
AS  THE  REQUIRED  RELIABILITY 
IS  THIS  VALUE  ACCEPTABLE?  (Y/N) 


PLEASE  PROVIDE  A  STARTING  BRACKET 
[B,C]  SEPARATED  BY  A  COMMA 
WHICH  YOU  BELIEVE  CONTAINS  A  ROOT 
THANK  YOU,  YOUR  STARTING  BRACKET  IS 


.iVv'V. 


0,1000 
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Computer  Output 

THE  OPTIMAL  VALUE  OF  *2  IS  =  250.421517 


THE  ORIGINAL  PROBLEM  WAS 

MINIMIZE  COST  5  *  N(  1  )+7*N(2)+6*N(3) 
S.T.  (1  -  ( . 4~N( !>)*(!  -  (.6~N(2))*(1  -  ( . 5^N( 3 > )  = 


THE  GP  SOLUTION  FOR  THE  PROBLEM  HAS  VALUES  OF 
N ( 1 )  =  4. 199369 
N( 2)  =  5.792162 
N(3 )  =  4.903512 

WITH  AN  OPTIMAL  COST  OF  90.96304798 
AND  A  RELIABILITY  OF  .91409482752 


THE  OPTIMAL 
N(  1  )  =  5 
N(2)  =  6 
N(3)  =  5 


INTEGER  SOLUTION  HAS  VALUES  OF 


WITH  AN  OPTIMAL  COST  OF  97 
AND  A  RELIABILITY  OF  .896896312 


